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1. Introduction 

Ideally, quantum computations are performed via transforming pure states of a physical 
system called 'quantum memory' to other pure states; see e.g. [H [21 [3l HJ O [6] . In this 
context, memory states transform unitarily. In most systems pure states may degrade 
to mixed states. This phenomenon amounts to decoherence. A well-known kind of 
decoherence is caused by extraneous influences unrelated to memory operations [6l[7]. 

Recently, Wu [HI El [10] introduced spatial variables in quantum computations by 
viewing the quantum memory as a scatterer: incoming particles are scattered from the 
memory and change its content. In this setting, unitary transforms apply to the com- 
bined system of memory and particles. The memory states do not transform unitarily 
unless the incoming signal is 'admissible'. In one space dimension, single-frequency 
waves are admissible [9]. In practice, however, incoming signals are pulses of finite 
duration. Thus, their use leads to additional decoherence, which we term 'impurity'. 
This kind of decoherence is connected specifically to memory operations, as was first 
discussed in [9]. The impurity of a two-state memory was analyzed via the relativistic [9] 
and nonrelativistic [H] Schrodinger equations. In the memory is allowed to interact 
with incoming particles only at one point by use of the pseudo-potential derived in [8] . 

In the present paper we extend the nonrelativistic formulation of [11] to reasonably 
general interaction potentials in one space dimension. Our starting point is to model 
the interaction potential as an imaginary quaternion [121 [131 [IH [IS] • In this formalism, 
the impurity measure of [TT] is expressed naturally in terms of an appropriate norm 
that depends on quaternionic commutators; see proposition I of section 12.21 In this 
context, pure states correspond to null hyperbolic quaternions [16]. For point interac- 
tions, scattering from the memory amounts to appropriate rotations of quaternions in 
the frequency domain. This approach offers additional insight into properties of the 
impurity measure used in [TT], and is amenable to computations. 

The paper is organized as follows. In section [2] we formulate the problem of impurity 
for a two-state quantum memory as a scattering problem with two coupled channels 
and general interaction potential in one space dimension: in section 12.11 we formulate 
the equations of motion by treating a local interaction potential as a quaternion; in 
section 12.21 we express the time evolution of the impurity measure used in [TT] in terms 
of commutators of quaternions; and in section 12.31 we describe an extension of this 
formulation to nonlocal interaction potentials. In section [3] we describe the general 
solution by invoking discrete schemes for amplitudes of suitable Fourier transforms in 
time. In section [4] we revisit the case with point interactions by use of the present 
formalism: in section 14.11 we focus on even wavef unctions; and in section 14.21 we treat 
odd wavefunctions. In section [5] we summarize our results and discuss open problems. 
Throughout the analysis we apply units with h"^/ (2m) = 1 where m is the particle mass; 
and denote quaternions by boldface symbols distinct from vectors. 
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2. Formulation 



In this section we describe the equations of motion for a general interaction potential. 
Subsequently, we derive an explicit formula for the impurity used in [Tlj . 



2.1. Equations of motion 

For a two-state quantum memory the field of the particle- memory system is the 
2x1 (column) vector 



'^{x,t) 



i'i{x,t) 
i)2{x,t) 



OO < X, t < +00 



(1) 



where iljj{x,t) (j = 1, 2) are scalar, square integrable functions in 1-1-1 dimensions. The 
vector field ip solves the Schrodinger equation, 

idtip = -dlifj - iq{x)ilj , (2) 

where — zq(x) represents the interaction potential and q is a 2 x 2 skew-adjoint matrix; 
see dH]) below. 

The q of (E]) is written as 

q{x) = J29''{x)aa, a = 1, 2, 3 , (3) 

a 

which we call an 'imaginary quaternion' [12]; g"'{x) are given real functions, aa '■= ica, 
and (/i = 0, 1, 2, 3) are the Pauli matrices with the usual convention ctq = 1, 

'nil fn ^ 1 ri n ' 

(4) 



0-1 

So, (Ta have the following properties. 

t _ ^2 



'0 


1" 




'0 -i' 




'1 " 


1 







_i 


, 0-3 = 


-1_ 



^eabc^c] aj^b, a, 6, c = 1,2, 3 



(5) 



where denotes the conjugate transpose (Hermitian conjugate) of a, and e^bc is the 
Levi-Civita symbol. Hence, we ensure that 

q \x) = — q(x) . (6) 

A simplified form of ([2]) is 

idt4> = -dli) - i[ge{x)\ie + go{x)\io]iJ , (7) 

where Ue^o are unit-length imaginary quaternions and ge,o{x) are even and odd functions, 
respectively. 

We note in passing that the 'discrete version' of ([2]) reads 



idti^ = —d1.ll) — ^^2qj(5(x — Xj)ip 



(8) 



where {q^ = q(a:j)} (j: integer) is a sequence of imaginary quaternions. By convolution 
of (IHl) with an appropriate kernel and for sufficiently dense partition {xj}, the resulting 
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solution can be arbitrarily close to the solution of ([2]). Because ([H]) is amenable to 
numerical computations, we discuss the relevant solutions in detail in section [31 

For later convenience, we introduce the Fourier transform in time of ip{x^t) by 
assuming that this signal contains only positive frequencies. With the definition 



oo 



Hx,t)= / — e-^-*0(x,u;), (9) 
Ztt 





the equation of motion ([2]) transforms to 

u(l){x,uj) = —dl(f){x,u) — iq{x)(t>{x,Lj) . (10) 

2.2. Impurity 

Next, we analyze impurity on the basis of ([2]). The reduced density matrix p{t) for the 
memory is the 2x2 matrix [TT] 

+ 00 

Pit) := / dx tP{x,t)iP\x,t) , (11) 



Q{x,t) := iptp^ 



which is obtained by tracing out the spatial variables. Of particular interest is the limit 
M := lim p{t) , (12) 

t— >+oo 

which is connected to the final memory state. For vanishing impurity, = M and 
tr(M^) = 1. In [11], the impurity measure is defined by 

Imp(M) := [1 - tr(M^)]^/^ . (13) 

In this subsection, we describe the M of (|T2l) . For this purpose, we form the density 
matrix 

Y^r'^Oa) =:-ir{x,t), (14) 

a 

where ip* is the complex conjugate of ipj, and the coefficients (// = 0, 1, 2, 3) are 

r' = ^xri^2 - , = - m') ■ (is) 

It follows that 

+00 +00 

p{t) = —i dxr{x,t) =: —im{t) , m(t) := / dxr{x,t) , (16) 



M = — ? lim m(t) =: —iuiout ■ (17) 
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An important remark is in order. With r = ir^crg + Xla^"^" (|T^ . the dual of r 
is defined by 

r'l-zrVo- J^r'^a, . (18) 



By virtue of flTSl) . we have 

rr^=:(r,r)=0, (19) 

i.e., r is identified with a null hyperbolic quaternion |16]; (■ , ■) is the Minkowski inner 
product by identifying r with the four- vector (r° , r^, r^, r^). Note that the m in ( IT6|) is 
an integral of null quaternions. Because of the convexity of the characteristic cone |14j . 
this integral produces either a time-like (mm'^ < 0) or a null (mm'^ = 0) quaternion. 
This property ensures that the 1 — tr(M^) in ( |T3l) is non-negative; see appendix. 

We now derive an equation of motion for Q{x,t). With the system Hamiltonian 
H := —Idl — ^q(a;), ([2]) becomes idtip = Hip. By (fT4l) . we readily obtain 

dtQ = -i{'H^if)if ^ + iifiifSH) , (20) 

where '?/'J^H denotes the action ofH onip ' from the left. Equation (|20|) is recast to the 
conservation law 

dtQ + d,P = [Q{x,t),q{x,t)], (21) 

where [Q , q] := Qq — qQ is the commutator of two quaternions and P is the 'flux 
matrix' 

1 



P = - 

i 



(dj)^^ -i^{d,^^) . (22) 



By (|T6l) and (|T71) . we need to integrate (1201) over space and then time. Accordingly, 
we obtain 

+ OD +00 

Am := uiout - ^in = J j dtdx [r(x,t) , q(a;)] , (23) 



— oo — oo 



where mj„ := limj^_oo m(t) is given. Consequently, 

M = -i(mi„ + Am) . (24) 

It is convenient to rewrite the Am of fl23l) in terms of the appropriate Fourier 
transform. Plancherel's formula [17] and definition f|T^ give 

+00 +00 

-^v{x,u), (25) 
Ztx 

-oo 

where r(x, uj) is the quaternion corresponding to the frequency-domain density matrix 
4>{x,u!)(p ^{x,uj) by ([9]), i.e., 

r{x,u) := i(j){x,iu)(f) \x,iu) , (26) 
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assuming that the signals have only positive- frequency content. Thus, we have the 
formula 

+00 +00 

— / dx [r(x, uj) , q(x)] . (27) 
27r J 



We note in passing that fl27j) can be formally generalized for higher space dimensions, 
where (j){x, uj) satisfies the vector Helmholtz equation. 

We now restrict attention to pure initial memory states. The incoming vector field 
has the product form 

V'm ~ tpinix,t)Sin t -OO , (28) 

where ipin{x, t) is the incoming particle wavefunction and Sj„ is the initial memory state. 
For example, we have [Uj 

did _,„,.,■. /777UI ^ 1, even wave 



V^.n(x,t)= / *-v^N/(o;)-<' ^' , (29) 

J ZTT sg(x), odd wave 

^ 

where sg(x) is the usual sign function, i.e., sg(x) = 1 if a; > 0, sg(a;) = — 1 if a; < and 
sg(0) = 0. The incoming quaternion is [llj 

+00 

£ \f{u)\'V^{s.r.S,l) . (30) 



We view ( 130|) as the general definition of f{uj) (without specifying any symmetry in x). 
By analogy with (fT9l) . mj„ satisfies 

mir^ini = . (31) 

Similarly, the condition tr(M^) = 1 (pure final memory state) corresponds to 
lUout^iut = 0- So, pure states are described by null hyperbolic quaternions. This 
mapping is one to one, as stated in proposition I below. 

The frequency profile |/(ci;)p in ( |30ll is chosen so that uiout is as close to null as 
possible. We state the following proposition. 



Proposition I. The impurity measure U^) reads 



Imp(M) = V2Km,„t, inout)\ = V2|2(mi„ , Am) + (Am, Am)| , (32) 

where {■ , ■) is the Minkowski inner product defined in / flgj) . (Thus, minimizing impurity 
is equivalent to minimizing the Minkowski norm corresponding to niout-) In addition, 
zero impurity is equivalent to niout being a null hyperbolic quaternion. 

A proof of fl32|) follows directly from the definition of v^out = "^M and f|T3|) . 
We sketch the main steps here; for details see appendix. By the representation of 
iWout in terms of the Pauli matrices we find | {mout , T^out) \ = detM. Recall that the 
impurity measure is Imp(M) = ^2 (detM) [llj, and for a pure initial state we have 
(mi„ , mi„) = 0. 
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Proposition I shows that definition (1131) for Imp(M) is a natural choice: the 
deviation from a pure state is expressed in terms of the 'length' of a hyperbolic 
quaternion. The impurity can be obtained from ( 132|) combined with (127|) . (l26l) and ( l30l) 
once is known. The reader is referred to section [3] for details on (j){x,uj). 

2.3. Extension 

In the case with a nonlocal interaction potential [8], the equation of motion is 

+00 

idtiJ = -dlif - / iq{x,y)ilj{y,t) . (33) 



The corresponding equation for Q = ipip ' is 



+00 



dtQ+d,P= / dy \^{x,t)^\y,t)q{x,y)~q{x,y)iJ{y,t)iJ\x,t)} .{34) 



If (i{x,y) = qV^(x)Va(y) and are scalar functions [8J, integration of fl34|) yields 
dtui = [vv , q] , (35) 
where the quaternion ry is defined by 

+0O +00 

rv= [ dyVp{y)^{y,t) I dx V^{x)i^ \x,t) . (36) 



3. General solution scheme 

In this section we describe the Fourier transform 0(x, tu) of %l){x,t) by solving (fTOl) . We 
start with ([H]), the discrete analogue of ([2]). The Fourier decomposition ([9]) reduces (IHl) 
to the form 

u;0(x, uj) = —d1(j){uj, x) — i5{x — Xj)qj , (37) 

j 

where 0^ := (f){xj,uj). The general solution of fl37|) is 

0(x,o;) = a+(cu)e*^ ^ + a_(o^)e-*v^ ^ + ^zqj0jG(x - x^) , (38) 

j 

where a± are reasonably arbitrary vectors and G{x, u) is the Green's function defined 
by 

G(x, to) = ( e*^l"l - e-^v^l^O . (39) 

By inspection of (1551) and fl5^ . (p reads 

0(x,cu) = cr+(x,cj) e*^^ + cr_(a;,cu) e"^^^ , (40) 
where only the values d±{xj, u) matter. 
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In principle, the coefficients d± can be determined from the incoming wavefunction 
ipin, formula fl28|l . In the limit t ^ — oo we have [11] 



do; 
2^ 



t — *• — oo 



(41) 



where the upper (lower) sign is taken for a; < (s > 0). In particular, we set 
X =Foo [Hj. So, the last relation gives c?+(x, ■) if x < — M+ and (i_(x, ■) if x > M_ for 
sufficiently large M± by use of the Fourier transform of ipin-, e.g. (l29ll . In the following, 
we find 0(x, •) everywhere via a scheme that determines d±{xj^uj). 

Next, we derive equations for d±j := d±{xj,uj) on the basis of (l37l) . For this 
purpose, we introduce the four-component vectors 

'd. 



d{x, Lu) 
Equation (1381) reads 



d. 



aiuj] 



^3 = ^-^ sg(j - A;) Vtkdi, 



where dj := d{xj,u) and Qj is the 4x4 matrix 



Wo- 



By taking differences in (H3!) we find the equation 

(/ + ^•+i)4+i = (^ - ^O^- , 
where / is the 4x4 unit matrix. By the identity (fij)^ = we write fH5l) as 

(^.^1 = (/ + ^],+l)-l (/ - vij)d,j = {i - n^+i) (/ - nj)dj . 

Let us assume that q^- have finite range, i.e. 
q,- = \j\>N, 
for some fixed positive integer A^. Define 

die ft ■= dj j < -N , dright ■■= dj j > N , 

which are constants. By introducing the 4x4 matrices 

i-i 

n 



N 



R, := (/ 



^0 



2Qu 



-N 



-N 



we derive the relations 



< N 



dright 



'] — -1 ^j "'left 

By the summation form (143|) we obtain the formulas 

die ft = a + ^jRj die ft, Rdleft = 



Rd 



left ■ 



y^^^ljRjd, 



left , 



(42) 
(43) 

(44) 

(45) 
(46) 
(47) 
(48) 

(49) 

(50) 
(51) 
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by which we find the relations 

I - R = 2'^njR 



d= -{I + R)dieft 



(52) 



(53) 



Thus, by (!50|) . diejt alone suffices to yield 0(x, u) in (!38|) . With recourse to (HTI) . the 
incoming wavefunction ipin furnishes immediately the 2x2 vectors d+^ieft and d-^right- 
By writing 

Ri ^2 

R3 R4 

where are 2x2 matrices, and using (!50|) we find 

d-,left = R4 ^{d^ .right — ^zd+,left) ■ (54) 

For the sake of simplicity, we assume that R4 is non-singular. The last relation completes 
the calculation of the four- vector die ft- Thus, we arrive at the following statement. 

Proposition II. Equation [3l\ ), with = for \ j\ > N, is solved by (E^j where a given 

by m, 

= d+j e'^ + d.j e-'^ , (55) 
andd±j are given by ^5^); d+^i^ft is determined by ( fTi] ) andd^^ieft is determined by (54^ - 



The continuous analog of (H3l) . which pertains to the solution of (ITOl) . is the 
Fredholm-type integral equation 



+00 



d{x,uj) = d{uj) - / dy sg{x - y)n{y,uj)) d{y,uj) 



where the 4x4 matrix Q is 

1 



f2(x, u) 



q(a;j 



w[x)q^[x) 
Ay/uj [—w*{x)q{x) — q(a;) 



w[xi 



(56) 



■ (57) 



(5^ 



Differentiation of (156!) leads to the Dirac-type equation 
dxd{x,uj) = —2Q{x,uj)d{x,u) . 

Note that the scheme underlying proposition II corresponds to solving ( |56l) by 
iterations. Because of the obvious connection of this scheme to the standard theory of 
integral equations [H], we do not discuss ( l56i) any further in this paper. Once (j){x, u) is 
known, Imp(M) can be calculated via proposition I in section [2^ In the next section, 
we apply proposition II to a delta-function potential [HI flT] . 

The procedure of this section, which applies to the Schrodinger equation ([2]) with 
a local interaction, can be extended to nonlocal interactions, equation fl33|) . but the 
algebra is more elaborate. In the next section we consider the case where the kernel of 
the interaction becomes a suitable pseudo-potential [Hf [TT]. 
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4. Point interaction 

In this section we revisit the theory of in tlie context of the present formalism, 
particularly of propositions I, II in sections 12.21 [31 For point interactions and a pure 
initial state of the memory, we calculate the impurity measure and show that the 
scattering amounts to rotations of quaternions in the frequency domain. We describe 
how a class of incoming finite-energy pulses can produce small impurity. 



4.1. Even wavefunctions 

First, we consider the Schrodinger equation 

idtip{x^t) = —dlip — iq6{x)ijj{0,t) , 
by which the particle interacts with the memory at the origin. We set 

q = ||u|| = 1 , 



(59) 



(60) 



i.e., u is a unit imaginary quaternion; ||u|| := u). The incoming wavefunction is 

assumed to be the even part of fl2^ [TT] . 

By virtue of (1271) and (l32l) . the impurity is measured in terms of the quaternionic 
commutator 

+ 00 



Am 



duj 
2^ 



(61) 



where r(0,a;) is 

r{0, to) = i(p{0, uj)$^{0, to) . (62) 

We now apply the formalism of section 131 in particular proposition II. In the present 
situation we have N = 0; the associated vector coefficients are d_i, do and di. By 
there is only one propagation matrix, i.e.. 



j^U 1 + j^U 



Thus, 



dleft — d-i — (/ + ^7o)'^0 , dright 

The introduction of the 2x1 vector ds by 
do = 



di = {I - Qo)do 



converts (EH) to 



di 



sg(j) 



9 



u 



-1,0,1. 



Vectors relevant to the Am of (16T1) are 

9 



di, 



d. 



,left 



1 + 



2v^ 



u]d^ 



(p{0,to) = 2ds 



(63) 
(64) 
(65) 

(66) 
(67) 
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By (1611) . the entanglement quaternion Am reads 

+00 +00 

Am = 4 / ^ [zdJJ ,gu]= [ ^ ^ e-'^[zd,J,i , u]e^" (68) 
J Zn J Zn 



where we conveniently defined the quaternion 

p(u;):=l + -^u=:pe^V (69) 



We apply the convention that the magnitude and phase are 



p{uj) = ^l + ^, ^(^)=arctan(^) . (70) 

We now simplify fl68|) by observing that dtu/dO = —{A/g)p'^uj'^^'^ and 

'idindii = \ f{^)?Un , ti„ := i{sinS^l) , (71) 
where Sj„ is introduced in f l29l) . Furthermore, we apply the identity 

A(e-ne^") = e-^"[t,u]e^", (72) 

where the operation e~^"te^" is a rotation which leaves the plane spanned by {(Tq? u} 
invariant. We can find two imaginary quaternions {v, w} that are orthogonal to u; 
then, let be the projection onto the space spanned by {v,w}. Equations (l68il - (172ll 
entail 

+00 

Am=-16y'^^^/VHP^{e-''"(PuUe^"} • (73) 


Integration by parts yields 

+00 

Am = 16 / ^ e-^(-)"(P.t„)e^(-)"(.;3/V(^)r)' , (74) 
J ^vr 


where the prime denotes differentiation with respect to uj. 

Example of small-impurity pulse. We proceed to describe how a class of incoming 
finite-energy pulses can produce an arbitrarily small impurity. Such pulses have of course 
a narrow spectrum sufficiently localized at a single frequency. 
A simple case of pulses with finite energy is described by 

{2TxY'^u^''^\f{uj)\'^ = H{uj -uq) - H{uj - Kuq) , fsT > 1, cjq > , (75) 

where H is the Heavyside function {H'{u) = and {ujo,K) are given parameters. 

The amplitude |/(a;)p is scaled by 4 J^"^ duj \f\'^y/uj so that the total pulse energy is 
fixed to unity. This normalization will be carried out in the impurity measure Imp(M) 
below. In view of (!74l) . we compute 

Am = 16(e-^^"P„t,„e^^" - e-^^"Put.„e^^") , (76) 
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where 

01 = arctanf-^) , 62 = arctanf ^ ; < ^2 < < ^ • (77) 

Note that the operation e^^^PutmC^" rotates Pu^m by —Id. 
With regard to mj„, by (!30|) we compute 



min = ^ij ^|/(^)Pv^ t,„ = 41n(K)t,, . (7^ 







We notice that Am hes in the plane spanned by {v,w} which is orthogonal to {l,u}. 
Substituting in (1321) and normalizing properly we find that the impurity measure equals 



Imp(M) = 4 / - cos(2gi - 262)} + ln(ir){cos(2g0 - cos(2g2)}| .^g. 

The right-hand side of this expression vanishes only for K = 1, but becomes arbitrarily 
small if — 1 ^ 1. In this limit, 

Imp(M) ^ ./I (Aoo) , A^:=|/^-l|.;o. (80) 

V 3 tuo + 5^ /4 

The behavior Imp(M) = 0{Auj) as Act; ^ is expected to be generic for any 
incoming pulse wavefunction that has spectrum sufficiently localized at ujo with support 
(bandwidth) of size Auj. The precise prefactor that enters the formula for Imp(M) 
depends on the specifics of the pulse spectrum. 

4-2. Odd wavef unctions 

Next, we turn our attention to the equation [H [11] 



idti^ = -dl^ + \J dy 6'^{y)4j{y, t) j , (81) 

where Sp{x) denotes ^'{x) modified to remove any discontinuity at x = from the 
function on which it acts [8]: 6'p{x)g{x) := 5'(x)[l — \imx^Q+]g{x) for a; > and 
6'p{x)g{x) := 5'(x)[l — lima^.^o-]5'(^) for x < 0. For simplicity we write 



hoo 



^,(0,t)= / dy6'Jy)^{y,t) . (82) 



—00 



The incoming wavefunction is assumed to be the odd part in (12^ [TT]. The Fourier 
transform in t of flHTl) gives 

ij$= -dl$-iciS'p{x)MO,u;) . (83) 

By proposition I in section 12.21 and equation fl23|) . the entanglement quaternion 
reads 

+00 

Am = j dt [n(0, t) , q] , (84) 
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where 

n(0,t) ■.= iMO,t)^JiO,t) . 
By use of the Fourier transform of t), we have 



5) 



+00 



Am 



da; 
2^ 



n(0, , q] , n(0, u) = i 4(0, uj)(j)J{0, uj) . 



The solution of (!83|) reads 
where Gx is the derivative of the Green's function (|39l) . i.e. 



1 
4 



G^fx, u) = sg(x) ( e^^^ + e-*^^ 



(86) 



(87) 



(88) 



Note that the pseudo-potential 6p gives zero when it acts on [8\. Thus, applying 6p 
on dHZD yields 

0,(0, 00) = ty^{a+{uj) - a-(cu))} . (89) 

Substitution of ( l89l) into ( IHTl) with = {d^, d-)^ leads to the formula 

l + fsg(x)v^u -fsg(x)v^u 



d{x, uj) 



According to the imposed antisymmetry, we introduce the 2x1 vector da in 













d- 




_-da_ 


have 










UJ) = 




d_{x, 




--{ 



(90) 



(91) 



1 - |sg(x)v/t^u[>4 . 



(92) 



Furthermore, by (|57|) . 

4(0,cu) =2«v^4 . (93) 

The vector din corresponding to in-states can be (i_(+oo, uj) or (i+(— 00, cj); compare 
to (HTj) . For example, 

4(^0^) = (l-|v^u)j; . (94) 

Hence, by analogy with section 14.11 it makes sense to define 

l + ^v^u=:p(a;)e''" , (95) 

where 

(96) 



viyj) = a/1 + ujg'^/4: , 6'(co') = arctan^ 2"^^ ' • 
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Equation (1861) for Am becomes 




(97) 





The assumption of a pure initial state amounts to using idind^^ from (17T]) . i. 
|/(co')ptj„ where tj„ = isinS^^ and Sj„ is introduced in fl28l) . Note that 
duo 4 ^ „ 





(98) 



Thus, by analogy with the symmetric case (section 14.11) we find the formula 



Am = 16 




{ 



e^(-)"(Put.„)e-^(-)"}(^=^/2|/(..)n' 



(99) 







The impurity Imp(M) follows by the procedure of section W7\\ 

A simple example of an incoming pulse wavefunction is described again by 
(27r)^^ci;^/^ \f{ijj)\'^ = H{ijj — ojq) — H(uj — Kcuq). The analysis for the impurity follows 
the steps of section 14.11 and is omitted here. Equation fl79l) should be recovered, where 
the angles 9i and 6*2 are now defined by 



5. Conclusion 

We introduced a general formulation of the nonrelativistic scattering from a two-state 
quantum memory in one space dimension. The key feature is to view the interaction 
potential as an imaginary quaternion. In the case with point interactions, scattering 
from the memory amounts to a rotation in the frequency domain of an appropriately 
defined incoming quaternionic state. 

We described the time evolution of the (entanglement) reduced density matrix in 
terms of the space integral of appropriate quaternionic commutators. By identifying 
quaternions with four- vectors we point out that, because of the space integration, the 
quaternions involved in the entanglement evolution are time-like. Accordingly, the 
impurity measure for the final memory state is described by a time-integral containing 
the Minkowski norm of time-like, hyperbolic quaternions. In the special case of narrow- 
band pulse wavefunctions, the resulting impurity Imp(M) is generically of the order of 
the pulse bandwidth. 

This work can be useful for addressing several questions. It is tempting to study 
successive scatterings from a quantum memory modeled by point interactions. In this 
case incoming states may not be pure but incoming quaternions are successively rotated 
in an appropriate sense in the frequency domain. An interesting question is how the 
impurity changes by this process. It is expected that Imp(M) always increases in such a 
case, especially if Imp(M) is thought of as 'entropy' in a sense [I9]. There is no rigorous 
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justification of tliis claim at tlie moment. Tlie connection of Imp(M) to interference 
effects critical to quantum computing sucfi as tfiose discussed in [20] was not addressed 
by our analysis. 

Another possible extension is the case of a n-state memory. A related issue is to 
define the appropriate algebra of n x n matrices that describe scattering in this context. 

Finally, it is interesting to consider relativistic massive particles within the present 
framework. A starting point would be the case of particles with spin 1/2. A perhaps 
naive model problem is the one- dimensional scattering from a two-state memory in the 
setting of Dirac's equation. In this case, the particle-memory system is described by a 
8x1 vector field. The study of this process by use of an analogous formalism is the 
subject of future work. 
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Appendix A. Proof of proposition I 

In this appendix we prove proposition I of section 12. 2[ In particular, we show that 
an arbitrary null quaternion can be written as a tensor product of the form ipip ^ . 
Furthermore, we show that the impurity measure Imp(M) is given by (!32|) . 
Suppose we have a null quaternion 

3 

r = irVo + 5^r'^<Ta , (A.l) 

a=l 

and let us write r := (r^,r^,r^) for the space vector. The 'nullity' property means that 
(^r°)^ = normalize so that r° = 1/2. The nullity is invariant under scaling. So, 

if we set h = 2r we have that h is also null with = 1 and \h\ = 1. An appropriate 
stereographic projection can identify the unit vector h with the complex number 

. :^ . (A.2) 

The inversion of this mapping yields 

1 + zz* ' 
{z - z*)/2t 

1 + zz* ' 

= . (A.3) 

1 + zz* 



r 



The substitution z := Zi/ Z2 in flA.Sp where -|- |2;2p = 1 gives 
r° = {zizl + Z2zl)l2 , 
= {zxzl + z\z2)l2 , 
= {z\z2 - zizl)/2i , 
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r' = {z,zl-Z2z;)/2. (A.4) 

These relations show that an arbitrary null quaternion can be written as the tensor 
product iptp ^ where 

2i|' + k2|' = l. (A.5) 

Next, we show formula fl32|) for Imp(M). By setting M = —ir, it is straightforward 
to calculate 

r2 = ( - (r0)2 - \f]-')ao + 2zr° J] r'^a, . (A.6) 

a 

By invoking the algebra of the Pauli matrices, we find 

1 - tr(M2) = 2[(r°)2 - \f\^] = -2(rr'^) , (A.7) 

so that 

Imp(M) = v^v^^ = V2y/\{r,T)\ = V2\\r\\ , (A.8) 

which confirms (l32l) of proposition I in section [2^21 In the above, r = {r^ , , r"^ , r^) is a 
four- component vector and < ■ , ■ > is the Minkowski inner product. 



Z2 
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